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Abstract 

We propose a cold atom implementation to attain the continuum limit of (1 + l)-d CP(A-1) 
quantum field theories. These theories share important features with (3 + l)-d QCD, such as 
asymptotic freedom and 0-vacua. Moreover, their continuum limit can be accessed via the mech¬ 
anism of dimensional reduction. In our scheme, the CP(A-1) degrees of freedom emerge at low 
energies from a ladder system of SU(A) quantum spins, where the N spin states are embodied by 
the nuclear Zeeman states of alkaline-earth atoms, trapped in an optical lattice. Based on Monte 
Carlo results, we establish that the continuum limit can be demonstrated by an atomic quantum 
simulation by employing the feature of asymptotic freedom. We discuss a protocol for the adi¬ 
abatic preparation of the ground state of the system, the real-time evolution of a false 0-vacuum 
state after a quench, and we propose experiments to unravel the phase diagram at non-zero den¬ 
sity. 
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1. Introduction 

Recently, there has been growing interest in developing physical platforms for quantum simu¬ 
lation of Abelian and non-Abelian gauge theories l[Tl|2]|3]|4][2|6ll7][3|9l. This effort is motivated 
by applications in particle and condensed matter physics, with the hope of developing quantum 
simulation Qol as a new tool to access regimes and phenomena complementary to, and beyond, 
classical simulations El. Previous work has focused on implementing quantum simulation of 
lattice gauge theories. An outstanding example is provided by cold atoms in optical lattices as 
a natural and controlled environment na, where the lattice gauge theory of interest emerges as 
a low-energy effective description of tailored atomic Hubbard dynamics || 2 ] 13 . Applications in 
particle physics, however, ultimately require taking the continuum limit, to eliminate artifacts 
due to space discretization. 
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Figure 1: a) Atomic setup for the implementation of the CP(A^-l) model, where SU(A^) spins are arranged on a 2-d 
bipartite spatial lattice of volume Lx L' with L:» L'. b) The N spin states are realized by AEAs occupying N < 21 + I 
hyperfine states, where I is the nuclear spin, with Zeeman splitting due to a uniform magnetic field, c) Identifying 
c' (cf. Eqs. § and l|^) on sublattice A, and «-* d on sublattice B, gives rise to the description of SU(A^) spins in 

the fundamental (red/dark), and anti-fundamental representation (blue/light), respectively, with 1 and N - I fermions per 
site, d) Interactions between SU(A^) spins T" and -7“" are generated via superexchange, where J oc f-jU . An energy 
offset y allows for full tunability of J. 


While some effective field theories emerge directly from cold atom systems in continuous 
space (i.e. without a lattice), here we construct lattice field theories from atoms in an optical lat¬ 
tice. Instead of following the standard procedure of Wilson’s lattice theory, where the continuum 
limit is approached by tuning a bare coupling constant IfTTII . we use the formalism of D-theory, 
in which the continuum limit emerges via dimensional reduction IfT^fT?! . 

We illustrate this idea for the relevant example of CP{N- 1) quantum field theories IIT4lfT5l . 
Such models have attracted interest in the context of particle physics as toy models for QCD, 
with which they share key features such as asymptotic freedom, the nonperturbative generation 
of a mass gap, and the existence of nontrivial 0-vacua inKia. In addition, in a condensed matter 
context CP(N- 1) models [Tfil have been discussed in relation to deconfined quantum criticality 

imiia. 

It has been shown that the (iH-l)-d CP(V-l) model emerges via dimensional reduction as the 
effective low-energy dynamics of certain (2-i-l)-d spin ladder models of S\J(N) quantum mag¬ 
netism im. In this paper we will show how this particular construction allows one to implement 
and approach the continuum limit of the CP(V-l) model in a natural and realistic way with 
fermionic Alkaline-Earth Atoms (AEAs) in an optical lattice Il20l 1271 l22l |2^ |2^ l25l l26l l27l 1^ 

|29l|3i 1211132 ■ 

The SU(A) spin model of interest, and the corresponding atomic setup are illustrated in Eig.[T] 
We assume fermionic AEAs with nuclear spin / representing SU(A) spins with N < 21+1. These 
atoms are loaded into a bipartite 2-d spatial lattice of volume L x L' (L '3> L'), realized as an 
optical superlattice. Such superlattices can be realized by superimposing a conventional optical 
lattice with lattice spacing a with a second lattice of lattice spacing V2a with a relative angle 
of 7r/4.The lattice depths and the interactions between the atoms are adjusted to achieve a filling 
with 1 and N - 1 atoms on the A and B sites of the bipartite lattice, respectively. Atoms on 
neighboring sites will then interact via superexchange (ie. Heisenberg-type) processes. 
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Our paper is organized as follows. In Sectionj^we will begin with a description of the spe¬ 
cific SU(A^) spin model of interest on a 2-d bipartite lattice. In Section]^ for self-consistency, 
we will review how the CP(A^-1) model emerges as the low-energy effective theory of the spec¬ 
ified SU(A^) spin model, following the discussion in Ref. El. We will then discuss how this 
construction allows for the continuum limit of the (l-i-l)-d CP(A^-l) model to be approached 
via the mechanism known as dimensional reduction. In Section|^we show how, by approaching 
the continuum limit in this way, we build a bridge which allows for their natural implementation 
via AEAs, while in Sectionj^we present a summary on the possible experimental imperfections. 
Finally, in Sections|^to|^ we describe the experimental signatures of the properties of the model, 
including asymptotic freedom and its relation to the continuum limit, the phase diagram of the 
system, and the quench dynamics of a false vacuum. 

2. Anti-ferromagnetic SU(A) Spin Model 

We begin by introducing the quantum spin model, from which the CP(A-l) model will be 
shown to emerge. As mentioned above, this model is defined on a 2-d bipartite lattice (see Fig.[^, 
and described by the Hamiltonian 

H^J Yj J>0, (1) 

(xy), xeA 


where and r” are generalized spin operators transforming under the fundamental and anti¬ 
fundamental representation of SU(A), residing on sites x e A and y e B of the even and odd 
sublattices, respectively. The spin operators satisfy [T", T*] = ^dxyfabcT^, where fate are the 
structure constants of SU(A), and the anti-fundamental representation satisfies T^ - -T^*. Note 
that in the case N - 2 the fundamental and anti-fundamental representations are equivalent, and 
Fq. Q reduces to the Heisenberg model. 

The system has a global SU(A) symmetry with total spin conservation. 


[H, T'] = \h, Y T'x 

xeA 



( 2 ) 


The Hamiltonian can also be formulated in terms of fermionic operators, by rewriting the 
spins as fermionic bilinears 


T" - V ja d 

^ X 2_J 

mm' 


- - S A"* d 

/ j ^xm^mm'^x, 

mm' 


(3) 


where A‘‘ are the generalized N x N Gell-Mann matrices, Tr[/l"/i^] = 2(5"*, and dxm annihilates 
a fermionic mode at position x, with the indices m,m' e {1,..., A) labelling N fermionic states. 
Beyond being of direct interest for CP(A-1) models, this class of Hamiltonians has been exten¬ 
sively discussed in the context of frustrated magnetism as a natural extension of the conventional 
SU(2) Heisenberg model (see, e.g.. Refs. caiia). In particular, in the 2-d limit, L' oo, it has 
been shown how the SU(A) symmetry undergoes spontaneous symmetry breaking ll^ for small 
N. We will exploit this property below while discussing the origin of the CP(A-l) quantum 
fields emerging from the spin Hamiltonian. 
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3. Connection to the CP(,N — 1) Model 


In this section we will review how the (l + l)-d CP(A^-l) model emerges as the effective 
low-energy dynamics of certain (2-i-l)-d SU(A^) spin models, showing how this theory undergoes 
dimensional reduction, approaching the continuum limit of the (l-i-l)-d CP(A^) model. These 
results have been reported earlier m. However, to keep the paper self-contained we review 
them here, before discussing how this leads to a natural implementation with AEAs. 

The CP(N- 1) models have been widely discussed in the context of low-dimensional quan¬ 
tum field theories. In contrast to their 0(N) counterparts, they display stable instanton solutions 
even for all N, as shown in Ref. m. This property, together with the fact that the models also 
show confinement and asymptotic freedom, are features that (l-i-l)-d CP(A-l) models share 
with(3-Hl)-d QCD. 

While the CP(A-l) model can be studied analytically in the large N limit lfT4l . here we 
show the emergence of the model for A = 3,4, allowing for quantum simulation complementary 
to analytical approaches. We begin with the SU(A) spin Hamiltonian given in Eq. Q. In the 
zero-temperature thermodynamic limit L, L' —» oo, the SU(A) symmetry breaks spontaneously 
down to U(A - 1) on, resulting in 2{N - 1) massless Nambu-Goldstone bosons described by 
fields in the coset space SU(A)/U(A - 1) = CP(A - 1). These fields can be described by A x A 
Hermitian projection matrices P, with Tr P - - P and P^ - P, and by the action 

S[P] = dx£ dyTr[p,d^Pd^P + ^d,Pd,P - ^{Pd,Pd,P - Pd,Pd,P)], (4) 

where is the spin stiffness parameter and c is the spinwave velocity. In the case of a finite 
extent L', as a result of the Mermin-Wagner theorem, massless Nambu-Goldstone bosons are 
forbidden; they pick up a mass m = 1/^, where ^ is the correlation length. Due to asymptotic 
freedom of the (1 -i- l)-d CP(A-l) model, ^ grows exponentially with L', 

^ oc expiAnL'pJicN)). (5) 

As L' increases, ^ becomes much larger than L', ^ » L', and the fields become independent of 
the spatial direction y. In this regime, the system thus undergoes dimensional reduction, where 
the dynamics can then be described by the effective action in the remaining space-time dimension 

mi, 

c 1 

5[P] = — J dt J dxTrld;,Pd;,P+^dtPdfPj-iffQ[P], (6) 

where Q{P} € n 2 [CP(A-l)] = Z is the topological charge, and - cf(L'ps) is the coupling 
constant of the dimensionally reduced theory. In terms of the lattice formulation, the vacuum 
angle is given hy 0 = nn mi, where n - L' la is the number of legs in the L' direction, and a is 
the lattice spacing. Eor N - 2 this reduces to the well-known 0(3) field theory description of the 
low-energy physics in Heisenberg antiferromagnets. 

We can now highlight an important consequence of this construction. In contrast to Wilson’s 
lattice field theory, in D-theory Il35l the continuum limit, ^ja oo, is approached by increasing 
L', not by decreasing a bare coupling constant. Due to the exponential dependence of ^ on 
L' = na, cf. Eq. the continuum limit is already approached for moderate values of L', which 
are accessible in current experiments. This strategy to regularize strongly coupled held theories 
is generally employed in the context of D-theory: in particular, in the D-theory regularization of 
QCD (3-i-l)-d gluon helds arise from dimensional reduction as collective excitations of (4-1-1 )-d 
SU(3) quantum links, while chiral quarks arise as domain wall fermions IH. 
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4. Implementation with AEAs 

We now show how the construction of the CP(A-l) model presented in the previous sec¬ 
tion allows us to implement and observe properties of a quantum held theory approaching the 
continuum limit in an optical lattice. 

The Hamiltonian Q is realized in a natural way in a system of AEAs trapped in an optical 
lattice, based on their inherent SU(A) symmetry m- Our implementation conceptually rests on 
two main ideas: hrst, using the formulation of the Hamiltonian Q in terms of fermionic degrees 
of freedom, as shown in Eq. 0. and second, implementing a particle-hole transformation to 
account for the fundamental/anti-fundamental representation of SU(A) spins with N > 3. In 
practice, this implementation exploits the toolbox already demonstrated in systems of trapped 
AEAs |l20l|2l]|22]|2l|23|25l|26l|27l|2a|29l|33[^ A number of experiments with Sr and Yt 
atoms have already been realized, where various strongly correlated phases have been observed, 
including SU(A < 6) Luttinger liquids ll30l . and SU(A) Mott insulators ||28]|36l. 

To begin, we consider a system of fermionic AEAs trapped in a 2-d bipartite optical lattice, 
cf. Eig. [T]a). We assume the 2/ H- 1 nuclear spin states to be split in energy due to a uniform 
magnetic held (Eig. [T]b)). The Hamiltonian of such a system is expressed in terms of localized 
Wannier functions ll^ as 


H = Hi + Hu, 

Hf — —t 

m (xy) 

Hu = Y ^Mx(nx - 1) + (7) 

X xeA 

Here Cxm is the annihilation operator for an atom with nuclear spin m 6 {-/,...,/) in the Wannier 
function localized at site x, and - 2,„ cl„,Cxm is the corresponding particle number operator. 
We denote by f the nearest-neighbor hopping amplitude, U is an on-site interaction energy, and V 
is an energy offset between the two sublattices. Note that here the scattering length is independent 
of the nuclear spin level m, providing the system of AEAs with a global SU(2/ -H 1) symmetry. 
The reason is that the electronic and nuclear spin degrees of freedom are decoupled, implying 
that two atoms in different Zeeman states will intereact equally, independent of the value of;«/. 
This emergent SU(A) symmetry has been shown to be accurate at a level of 10“* ll29l . We also 
note that one can take N < 2/ -H 1 by initializing atoms into a subset of the magnetic states; 
if some Zeeman states are initially empty, their initial populations will not change due to the 
absence of spin-changing collisions. AEAs can realize SU(A) physics up to A = 10 (e.g. ®^Sr 
and ^^^Yb), but here we will concentrate on the CP(2) model, i.e. A = 3, since it is particularly 
interesting from a theoretical viewpoint. 

We initially occupy each site of the A sublattice with 1, and each site of the B sublattice 
with A - 1 atoms. While the A sublattice spins in the fundamental representation of SU(A) 
are embodied by a single fermion, the B sublattice spins in the anti-fundamental representation 
are embodied by A - 1 fermions, which are equivalent to a single hole, cf Eig. [T]c). Using 
the fermionic representation of Eq. we identify Cx with dx on sublattice A, and with t/] on 
sublattice B. 

We now consider the Hamiltonian Q in the strong coupling regime, f <si U,V: the contri¬ 
bution of H, causes virtual tunneling processes within the subspace of states with hxed particle 
numbers per site (the eigenstates of Hu), thus generating SU(A) superexchange terms BTlI . In 
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Figure 2: Panel a) Energy of the ground state (red, solid) and excited state (blue, dashed), in units of the spin coupling J, 
as a function of the adiabatic parameter r. Schematics show the evolution of the system from t = 0 (top left) to r = 1 (top 
right). Panel b) Time evolution of the false vacuum state for a 14 site system. Upper panel: Order parameter calculated at 
each bond as a function of time. Lower Panel: Time evolution of the order pai'ameter summed over all bonds (see text). 


this regime, to second order in tjU, the Hamiltonian Q emerges with ifT^ 


(-y + u{N - 3))(y - u(N - 1)) ’ 

where an antiferromagnet requires J > 0. 

The ground state of the system can be prepared via an adiabatic protocol; we start by prepar¬ 
ing a band insulator with N particles per site on a simple square lattice. The population in each 
spin state can be controlled using, e.g. optical pumping ll^lSTIl . Each site is subsequently split 
into a double-well by adiabatically ramping up a superlattice, realizing a system of generalized 
SU(A^) singlets akin to what has already been realized using bosonic alkali atoms Eli. 

The barrier between the wells is then adiabatically turned off, realizing the full quantum 
dynamics of Eq. Q. This procedure relies entirely on existing techniques, and can be applied 
for = 3,4 and various L'. For a single chain this works as follows: starting from a perfectly 
dimerized initial state, the inter-well exchange is switched on according to the time-dependent 
Hamiltonian 

h{t) = -(1 - t)7 ^ + th, (9) 

xeA 

where re [0,1] is the adiabatic parameter. The corresponding low-lying spectrum of a 14 site 
system is shown in Fig.j^a). The system does not undergo a phase transition during this process; 
the gap does not close while changing t from 0 to 1. [^ This ensures that an adiabatic ramp can 
be performed on time scales shorter than 1/7. 

5. Experimental Imperfections 

In this subsection we analyse how robust our proposal is in light of the possible imperfections 
present in such systems of AEAs in an optical lattice. The main sources of imperfections are; 

• effects of atom losses. 


’ We have numerically checked that the gap remains of order 0.67 even at L = 60. 
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• imperfect loading of the initial Mott phase, 

• role of external confinement and control over the number of legs, 

• spin population imbalance. 

Below, we discuss in detail the relevance of each of these points for the observability of CP(A^-1) 
physics within our scheme. 

5.1. Atom loss 

One of the most prominent sources of imperfection in an optical lattice experiment is the loss 
of atoms from the lattice. Such losses can occur both while preparing the initial Mott insulator 
phase, and throughout the evolution, due to, for example, three-body collisions. 

If we focus on = 3, sites occupied with 3 atoms only occur virtually, and thus three-body 
collisions are not a relevant source of imperfections. 

While sites with 3 atoms do occur at the initial stage of our adiabatic state preparation, the 
lattice at this point has twice the lattice spacing, which results in much more localized Wannier 
functions. 0 

5.2. Imperfect loading of the Mott state 

We consider in more detail defects in the initial Mott insulator phase, which arise in the 
preparation of the band insulator with N atoms per site. Recent experiments have reported such 
errors on the level of 1 missing atom per about 100 sites for bosonic systems 1391 EOl. Similar 
numbers have been reported for SU(2) fermions ll4Tll42l . while for SU(A^) {N > 2), the additional 
effect of Pomeranchuk cooling would probably lead to an even smaller number of defects Il43l . 

In our setup, this error implies the following; When we split each initialized site to create the 
required double-well structure, we do not have 1 and - 1 atoms on each sublattice necessary 
to realize an SU(A^) spin in the fundamental and anti-fundamental representations, respectively. 

For concreteness we focus on A^ = 3 in the remaining discussion. In this case, the Mott insu¬ 
lator phase should be initialized with N - h atoms per site and the anti-ferromagnetic condition 
requires that V < 2U. If one atom is missing, for fixed U, there are two possible configurations 
when the site is subsequently split into a double-well, corresponding to V < U and V > t/. In 
the first case, where V < U, one atom will be present on each site. In terms of the spin descrip¬ 
tion of the system this imperfection corresponds to having a fundamental spin on a site where a 
spin in the anti-fundamental representation should be located. The spin-spin coupling between 
two spins of the same representation is 7 = fiUI(V^ - U^). In the second case, where V > U, 
the ground state will have both atoms on the same site, with the other site empty. In the spin 
description, this corresponds to one spin missing from the system. 

In general, the system’s dynamics in the presence of such imperfections would be described 
by a generalized SU(A^) t-J model. However, for the regime 2U > V > U (which is experi¬ 
mentally reachable and satisfies the anti-ferromagnetic coupling constraint 7 > 0 in Eq. (7)), the 
dynamics of the system simplifies. In this case, the impurities are ‘static’, as moving atoms in 
those partially filled sites is an off-resonant process. The dynamics of interest is then described 


^Assuming a loss rate of A's ^ 10 typical for fermionic isotopes of AEAs ED. an estimate of losses for a 

lattice spacing of ~ \fim and a ramp time of 20ms, is 0.2% This is in contrast to 13% for a lattice spacing at half of this 
(~ 500nm). 
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by the SU(A^) spin ladder model, with a variable percentage of defects, depending on the initial 
number of missing atoms. We remark that the effects of such ’non-magnetic’ impurities have 
been investigated in a related scenario for strongly frustrated magnets as well l44l . 

5.3. Role of external confinement and a fuzzy’ value of L' 

In our description of the system we have relied on two assumptions about its geometry, 
namely i) that the system is confined in a region of constant density, and ii) that L' = na was 
given by a well-defined value of n, the number of legs in the L' direction. The latter is crucial for 
a well-defined value of 0 in the action of the CP(A^-l) model, thus this point is essential for the 
proper implementation. 

In current experiments there are two ways to control ’sharp’ boundaries. The first is to employ 
a box potential on the optical lattice which allows for a sharp edge of the system. For system 
sizes ~ 5Qfj.m the boundary effects will affect atoms within ~ 2fj.m of the box boundary Il45ll46l . 
For Bose-Einstein condensates, the effects of such boundaries have been quantified to correspond 
to extremely flat potentials of the form V{r) oc ((r - rtfijaY with a > 10, and cq ^ 25/rm. The 
effect on the spin Hamiltonian is an adjusted coupling, J, between the spins, due to an effective 
chemical potential in the presence of the box potential — an extremely small effect in our case. 

The second alternative is provided by the rapidly expanding technology of quantum gas mi¬ 
croscopes (see, e.g.. Ref. iQi for a review). In this case, sharp boundaries in the system can be 
imprinted at the single-site level using blasting beams, so that the system dynamics is effectively 
confined into sharp boxes. The additional underlying confining potential will play no role in this 
case, as it does not affect the spin degrees of freedom. 

5.4. Spin population balance 

One further imperfection which could be present in the experimental realization is the ini¬ 
tialization of the system with unequal numbers of atoms in each of the N spin states. In reality, 
these spin states can be controlled via optical pumping, with current experiments achieving an 
accuracy of less than 1 error in ~ 100 sites BTI . In the case of an error, the low-energy effective 
model is still the CP(A^- 1) model, however, with a non-zero particle number (see the discussion 
in Sec.|^. We remark that using spin-dependent in situ imaging (as in Ref. 1481 1 could also be 
useful as an efficient post-selection method to detect possible defects. 

With this in mind, the possibility of selecting spin populations using optical pumping is 
actually available as an additional feature of our implementation: By intentionally loading an 
imbalance of spin states, we can investigate the CP(A^ - 1) model at non-zero chemical potential. 
One can also employ Monte Carlo methods to investigate this problem, which will be the subject 
of future work. 

6. Continuum Limit 

In order to demonstrate explicitly that the SU(3) spin ladder gives rise to the CP(2) model 
in the continuum limit, it is vital to study the correlation length and verify that it increases 
exponentially with the size L' of the extra dimension, cf. Eq. By means of Monte Carlo 
simulations with a loop cluster algorithm |491 we have calculated the spatial correlation length f. 
I^Eor even L' we obtain Eig.j^a), which indeed shows the anticipated exponential increase of f 

^The Monte Carlo simulations were performed with periodic boundary conditions in the L direction, and open bound- 
aiy conditions in the L' direction. We measured the correlation length f and the second moment correlation length ^2 to 
ensure that their difference is negligible. 
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Figure 3: Results for the correlation length ^ (in the L-direction) from Monte Carlo simulations at L = 1500^. Panel a): 
Results at fSJ = 10 (blue triangles), pJ = 100 (black circles) and= 1000 (red diamonds). When the correlation length 
becomes comparable to the relevant energy scale set by the temperature it begins to deviate from the form in Eq. 
Panel b): Results for the correlation length ^ at jSJ = 100 in the presence of 0% defects (blue triangles), 0.1% defects 
(red diamonds), 0.5% defects (black crosses) and 1% defects (green squares). 


with L', in agreement with asymptotic freedom of the CP(2) model that emerges via dimensional 
reduction. The increase of the correlation length for L'/a = 4 to 12 should already be accessible 
using current experimental techniques. Even at temperatures around /3J ~ 10, corresponding 
to a physical temperature of the order of nK, the correlation length f(6a) is close to 10a, and 
already falls on the exponential that indicates asymptotic freedom. The correlation length ^ can 
be measured in a cold atoms setup via Bragg spectroscopy or through noise correlations Col. 

6.1. Continuum Limit in the presence of imperfections 

After the previous qualititative discussion on the effect of imperfections in the implementa¬ 
tion with AEAs, here we consider quantitatively the effect of such imperfections on the observa¬ 
tion of the correlation length f. The key question we would like to answer here is, whether or not 
asymptotic freedom persists in the presence of such defects, and at which point the exponential 
growth of f sets in. At a qualitative level, one would expect that the phenomenon stays intact for 
low concentrations, since a small number of defects will not drastically affect the (2-i-l)-d limit, 
where spontaneous symmetry breaking occurs. Nevertheless, numerical simulations have been 
performed in order to make concrete statements regarding experimental realizations. We have 
performed a Monte Carlo study in which defects have been modelled as an empty lattice site 
and 50 realizations of randomly distributed defects have been averaged over. We investigated 
concentrations of 0.1%, 0.5% and 1% for j3J - 100. In each case there is an L' range in which 
the exponential law demonstrating asymptotic freedom is clearly visible. 

One can see from the results in Eig. 3.b) that the correlation length at a given L' increases with 
the defect concentration, thus the regime where dimensional reduction occurs sets in earlier. As 
well, independent of the defect concentration there is a scale — set by the inverse temperature — 
at which the correlation length saturates due to thermal fluctuations. This behavior is noticable 
in Eig. 3 both in the case with and without defects. Because the correlation length at a given L' 
increases with the defect concentration, this implies that the L' range in which an exponential 
behaviour is visible shrinks with the defect concentration. We have conhrmed that the large 
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L' behaviour of the curves in Fig 3.b) is a thermal effect by performing some simulations at 

PJ = 1000 . 

Overall, the only quantitative change defects cause is to renormalize the spin stiffness of the 
system such that it increases with the defect concentration. This is illustrated by the exponential 
fit lines in Fig 3.b), where one can see the gradient of these lines increasing with defect concen¬ 
tration. The data are fitted up to values of L' where this temperature effect becomes prominent 
which, for pj - 100, is on the order of ^ ~ 20a. 

1. Finite Density Phase Diagram 

Just like QCD, CP(A^-1) models have a finite density phase diagram that is worth exploring. 
While in QCD a chemical potential /i can be coupled to the baryon number, in the CP(2) model 
two chemical potentials, /rs and can be coupled to the global SU(3) symmetry, which is 
thereby explicitly broken down to U(l)xU(l), or to SU(2)xU(l) along the solid lines in Fig.j^a). 
It is then interesting to ask whether the U(l) symmetries are affected by Berezinskii-Kosterlitz- 
Thouless (BKT) transitions. In a cold atom experiment a finite density situation corresponds 
to loading the optical lattice with unequal numbers of atoms in the three Zeeman states, cf. 
Subsection |5.4[ corresponding to additional terms in the Hamiltonian, 

- Hid\A^d^, ( 10 ) 

where and are the diagonal Gell-Mann matrices. A BKT transition is then signalled by the 
‘condensation’ of bosonic molecules, formed from two fermions (thus forming a bosonic pair) 
in specific combinations of Zeeman states, in the spirit of color superfluidity ifSOllSTlIS^ . 

8. Spontaneous C-breaking iA6 = n 

Having advocated the feasibility of approaching the continuum limit in a quantum simulation, 
we now consider an odd number n of transversely coupled chains, corresponding io 9 - ji. At 
this point, analytical considerations suggest a first order phase transition with spontaneous charge 
conjugation (C) symmetry breaking Il53l . which has been confirmed numerically m- In our 
proposed experimental realization with discrete spins, C corresponds to a shift by one lattice 
spacing in the longitudinal 1-direction, —> -T®*., -T^* where the sites x and x + 1 

belong to the A and B sublattice, respectively. An order parameter which signals C-breaking in 
spin systems is given by ll^ 


D = - r^rc), (ii) 

xeA 

which, equivalently, detects dimerization. When C is preserved (n even, 9 - 0) D vanishes, 
whereas when it is spontaneously broken (n odd, 9 - n) there are two degenerate ground states 
with opposite non-zero values of D, see Fig. In a cold atom setup, measuring the singlets 
contributing to D has been proposed and demonstrated via spin-changing collisions mm. A 
possible adiabatic preparation scheme for the realization of such generalized resonating valence 
bond states is illustrated in Fig.j^a). 


10 



Figure 4: Panel a) Conjectured phase diagram of the CP(2) model as a function of the chemical potentials /is,/is- In the 
vicinity of /i 3 = /ig = 0, the system is in the vacuum state. Besides a normal fluid phase (yellow/light), high density color 
superfluid phases (red/dark) are expected to appear. The phase diagram can be explored in the proposed cold atom setup. 
Panel b) Symmetry of the ground states for 6 = n. Above: In the ground state bonds emerge between two neighboring 
sites on sublattice A and B. The ordering of the bonds gives rise to a double degeneracy of the ground state. Below: 
Cartoon of the two degenerate ground states. A solid line indicates a dominant while no line represents a 

smaller value. 


9. Quench Dynamics Decay of a False Vacuum 

The possibility of initializing states composed of singlets provides an opportunity to investi¬ 
gate real-time quenched dynamics driven by H. In a finite system with open boundary conditions 
and even L, one can investigate how a false vacuum |-), cf. Fig.|^b), decays as a function of time 
after the Hamiltonian H is switched on. Such a false vacuum decay can mimic processes in in¬ 
flationary early universe cosmology, as well as bubble nucleation at a first order phase transition. 
Due to energy conservation, the false vacuum, which has an energy cost for any finite system size 
L, cannot decay fully into the true vacuum |-i-). Instead one expects damped coherent oscillations 
between the two vacuum states. In order to quantify the false vacuum decay in real time, for a 
single chain we consider the order parameter ll^ 

D(t) = - r^r-pi'PCO), (12) 

xeA 

which indicates whether singlet states predominantly form on the even or odd bonds. This order 
parameter is maximal for |-) (all even bonds have a singlet). In Fig.j^b), we show the dynamics 
of the false vacuum decay evaluated by exact diagonalization of a L/a = 14 site system starting 
in the initial state |-). At times f « 1/7 the even singlets (blue) are stable. At later times, the 
false vacuum decays, with correlations remaining only in the central part of the system, while the 
bonds close to the boundary revert the order. The decay of the full order parameter is depicted 
in the lower panel, which indeed shows coherent oscillations. Both the order parameter and 
the local singlet projectors can be experimentally measured as discussed in the previous section. 
While moderate system sizes can be reached using exact diagonalization, the real-time dynamics 
in the continuum limit is inaccessible to classical simulations. Experiments using the present 
scheme would shed light on the real-time dynamics of false vacua in CP(A-1) models. 
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10. Conclusions 


We have outlined a proposal for the quantum simulation of a CP(A^-1) quantum field theory 
using cold atoms trapped in an optical lattice with a ladder geometry. Our work shows how the 
continuum limit can be assessed using dimensional reduction, and how paradigmatic phenomena 
such as asymptotic freedom can be observed in cold atom experiments. Extending such inves¬ 
tigations to non-Abelian gauge theories would provide an indispensable tool for the quantum 
simulation of fundamental theories such as QCD at finite baryon density. 
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